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Abstract 

We present numerical evidence for the existence of spinning generalizations for non-topological 
Q-ball solitons in the theory of a complex scalar field with a non-renormalizable self-interaction. 
To the best of our knowledge, this provides the first explicit example of spinning solitons in 3 + 1 
dimensional Minkowski space. In addition, we find an infinite discrete family of radial excitations 
of non-rotating Q-balls, and construct also spinning Q-balls in 2 + 1 dimensions. 
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I. INTRODUCTION 



Solitons are important ingredients of models in high energy physics. Apart from be- 
ing responsible for various non-perturbative quantum phenomena, solitons are interesting 
in themselves, since they can be viewed as field theoretic realizations of elementary parti- 
cles. It is from this viewpoint that solitons were originally introduced into physics in the 
context of the Skyrme model - as models of hadrons. The subsequent developments have 
revealed soliton solutions in many other non-linear field theories in Minkowski space, such 
as monopoles, vortices, sphalerons, Q-balls, etc. These solutions describe localized, particle- 
like objects with finite energy. Their spectra of energy and charge are typically discrete. In 
addition, these solutions are regular everywhere, a property which is especially appealing. 
So far, however, solitons have been lacking one important feature of elementary particles: 
the intrinsic angular momentum J, which is zero for all known classical solutions. 

It is sensible to ask whether stationary rotating generalizations for the known static 
soliton solutions exist. More precisely, one is interested in finite energy, globally regular, 
non-radiating solutions for which the spatial integral of the T 0ip component of the energy- 
momentum tensor, 



is non-vanishing. In this definition the condition of stationarity (absence of radiation) is 
important. Indeed, it is always possible to construct a field configuration such that J ^ at 
the initial moment of time. Physically this would correspond to exciting the soliton to give 
it an angular momentum 1 . However, when the time evolution starts, the received excitation 
will be most probably immediately radiated away, leaving behind a non-rotating object. 

When talking about rotating solitons, it also seems sensible to distinguish between two 
types of rotation: spinning and orbiting. Spinning is associated with the intrinsic angular 
momentum, in analogy with the quantum-mechanical spin. Classical spin excitations, if 
they exist, should live in the one-soliton sector because they are excitations of an individual 

1 In the literature one can often find explicit examples of 'solitons in the rigid rotator approximation', for 
which the integral (|l|) does not vanish, as for instance 'rotating' Skyrmions (!]], knots etc. These 
configurations, however, are not solutions of the equations of motion, and at best they can be viewed as 
the initial values for the dynamical evolution problem. 




(1) 
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object. One expects that the corresponding angular momentum will assume only discrete 
values. Solutions describing spinning solitons in Minkowski space in 3 + 1 dimensions are 
not known. For example, it is not known whether the 't Hooft-Polyakov monopoles can be 
given classical spin 2 . In fact, it has been shown that, at the perturbative level, monopoles 
do not admit stationary rotational excitations in the one-monopole sector This means 
that monopoles cannot rotate slowly, with | J\ <C 1. Monopoles with finite (discrete) values 
of J are not yet excluded. However, to decide whether such solutions exist requires to go 
beyond perturbation theory and solving the complete coupled system of the Yang-Mills- 
Higgs (YMH) partial differential equations (PDEs), which is an exceedingly difficult task 3 . 

On the other hand, one can consider relative orbital motions in composite many-soliton 
systems. For example, one can imagine a rotating soliton-antisoliton pair balanced against 
mutual attraction by the centrifugal force. The spectrum of the angular momentum is 
then expected to be continuous. Solutions describing such orbiting solitons can actually be 



constructed. In the case of monopoles, for instance, there exist [23, 24] solitonic solutions 



of the YMH field equations with vanishing Higgs potential in the sector with zero monopole 



charge. Such static, purely magnetic solutions have been explicitly constructed in they 
describe monopole-antimonopole pairs balanced by a repulsive force of topological nature. 
Now, there is a simple way to add angular momentum to these solutions [ I I ] by using 



the global symmetry of the field equations which mixes the Higgs field H and the electric 
potential $: 

H — > H cosh 7 + $ sinh 7 , $ — > <D cosh 7 + H sinh 7 . (2) 
Applying this transformation with an arbitrary 7 7^ leads to solutions with an electric field. 



It is important to note that this also changes the angular momentum to J ~ sinh 7 [11 



1 



2 There is also the possibility to associate the spin of the monopole with that of fermionic zero modes living 
in the monopole background. In this case, however, the spin is not classical, and in fact is not related to 
the monopole itself. 

3 It was argued in |2f| that at least within the minimal axial ansatz, rotating 't Hooft-Polyakov monopoles 
can be excluded also at the non-perturbative level. 

4 This trick works only for non-BPS solutions - those obeying the second order YMH equations but not 
the first order Bogomol'nyi equations. For BPS solutions the angular momentum is invariant under (^|). 
For this reason one cannot use this method to produce spinning monopoles, because non-BPS solutions 
with finite energy and unit topological charge are not known |L9|. 
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(in fact J = Q, where Q is the electric charge [p5|j). These new solutions can be interpreted 
as describing the system of a monopole and antimonopole rotating around their common 
center of mass. 

Another example of systems which could be classified as orbiting solitons are rotating 
vortex loops. In certain models there exist vortices with stationary currents along them; 
'superconducting vortices'. One can argue || [7|] that taking a finite piece of such a vortex, 
bending and closing it to form a loop, finally adding a momentum along the loop, leads 
to an object (vorton) described by a stationary solution of the equations of motion 5 . The 
angular momentum in this case is associated with the macroscopic circulation of the current 
along the loop. This type of motion could be naturally classified as orbital rotation. 

On the other hand, one can also consider the rotation of a straight vortex along its 
symmetry axis. This would correspond more closely to the notion of an intrinsic spinning 
rotation. Although one can show that for the Nielsen-Olesen vortex such spinning excitations 



do not exist, they can actually exist in other models [g, [10], [12], [13], |2TJ. However, this only 
gives spinning solitons in 2 + 1 dimensions, while their 3 + 1 dimensional analogs will have 
infinite energy due to the infinite length of the vortex. 

In summary, we are not aware of any spinning solitons in Minkowski space in 3 + 1 
dimensions. At the same time, such solutions are known in curved space. Indeed, there are 
many rotating solutions in General Relativity. In the pure gravity case they comprise the 
family of Kerr-Newman black holes. These are very similar to solitons, but they are not 
globally regular and contain a curvature singularity. There are also globally regular rotating 
solutions p0| , but these require a source of a non- field theoretic origin. Interestingly, there 



exist two explicit examples of gravitating spinning solitons in pure field systems. 

The first example is provided by rotating boson stars |22], ^6|. These are solutions for 
a gravity-coupled massive complex scalar field with harmonic dependence on time and on 
the azimuthal angle, $(t,r, v?) = exp(iut + iN(p)f(r, with N integer. The energy- 
momentum tensor is time-independent, and the Einstein equations together with the Klein- 



Gordon equation for /(r, admit globally regular, stationary particle-like solutions |]22 
The angular momentum is quantized as J ~ N. The solutions with N > can be regarded as 
spinning excitations of the fundamental static, spherically symmetric solutions with N = 0. 



To our knowledge such solutions have not been constructed explicitly. 
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This example is instructive in the sense that it is clear 'what rotates': this is the phase 
ut + Nip. In this connection it is worth noting that the very notion of rotation in pure field 
systems is very different from that for ordinary rigid bodies. Indeed, it is meaningless to 
say that a given element of volume of a field system actually 'performs revolutions' around 
a given axis. Of course, one can imagine a solitonic object with a field perturbation running 
around it. However, such a moving perturbation will be most probably immediately radiated 
away. The example of boson stars thus shows that one can nevertheless have a rotating phase 
which is not radiated away. 

In stationary rotating systems without explicit time dependence the rotation will rather 
be associated with certain non-linear superpositions of the multipole moments of the fields. 
For example, in systems with vector fields, angular momentum may be present due to a 
non- vanishing integral involving the Poynting vector, 

J — J fx (Ex B)d 3 x . (3) 

A stationary, globally regular configuration for which this integral is non-zero would corre- 
spond to a rotating soliton. 

It is possible that such solutions could exist for the SU(2) Yang-Mills fields coupled to 
gravity. The Einstein- Yang-Mills (EYM) field equations admit globally regular, particle like 
solutions [Q. These gravitating EYM particles are static, spherically symmetric and neutral 
(their purely magnetic gauge field strength decays asymptotically as 1/r 3 ). Now, for these 
solutions one can perturbatively construct stationary, globally regular, axially-symmetric, 
slowly rotating generalizations Q . Surprisingly, their spectrum of J is continuous, as if they 
were composite objects, which is presumably due to the special feature of the field system 
consisting of only massless physical fields. Unfortunately, it is not clear at the moment 
whether these perturbative solutions exist also at the non-perturbative level - the analysis 



of [|T^, |25] indicates the opposite. It is however still possible that the solitons exist, but 
within a more general ansatz than that considered in |T5| , f25f . 

In summary, spinning solitons have been found only in curved space. It is then sensible 
to ask if spinning solitons without gravity exist at all. In principle, it is not excluded that 
only gravity can support the relevant rotational degrees of freedom. In order to rule out 
this logical possibility, we have undertaken an attempt to construct spinning solitons in flat 
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space. 



The solitons we have chosen to 'rotate' are Q-balls |5], [18]. These are solutions for a 
complex scalar field with a non-renormalizable self-interaction arising in some effective field 
theories. These solutions circumvent the standard Derrick-type argument due to having a 
time-dependent phase for the scalar field. In this sense they are somewhat similar to the 
boson stars. In the simplest spherically symmetric case the fundamental Q-ball solution 
was described by Coleman ||. Dynamical properties of these objects have been studied in 
@. Q-balls also appear in supersymmetric generalizations of the standard model, where 
one finds leptonic and baryonic Q-balls [[l(J which may be responsible for the generation of 



baryon number or may be regarded as candidates for dark matter [17 



In the present paper, we first recall the properties of the fundamental spherically sym- 
metric Q-balls. In addition, we present an infinite family of their radial excitations which 
have not been reported in the literature before. We then turn to rotating solutions and first 
consider them in 2 + 1 spacetime dimensions, where the problem reduces to an ordinary 
differential equation. Finally we consider the full problem in 3 + 1 dimensions and explicitly 
construct solutions with angular momentum. To our knowledge, our analysis gives the first 
explicit example of spinning solitons in flat space. 



II. THE MODEL 

Let us consider a theory of a complex scalar field in 3 + 1 dimensions defined by the 
Lagrangian 

C = d^d^* -U(\&\) . (4) 

It is assumed that the potential U has its global minimum at $ = 0, where U(0) = 0, while 
U — > oo for |$| — > oo. In addition, the potential must fulfill a particular inequality (Eq.(|T5|)) 
which will be discussed below. The potential may also have local minima at some finite |$|, 
as is shown in Fig. [I], but this is not necessary. 

The global symmetry of the Lagrangian under $ — > <&e ia gives rise to the conserved 
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FIG. 1: The qualitative shape of the poten- 
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FIG. 2: The effective potential V{4>) = 
±wV - U(4>) in Eq. dH). 



charge 

Q = iy d 3 a;($*$ - . (5) 

The fundamental Q-ball solutions of the theory are minima of the energy for a given Q 
0. Since $ should depend on time for Q to be non- vanishing, one assumes that $ has a 
harmonic time dependence. In the spherically symmetric case, 

$ = (f)(r) e luJt , (6) 

where <p(r) is real. The potential {7(|$|) = U(<p) and the energy- momentum tensor, 

Tp, = c\$* + d v § d^* -g^C, (7) 

(g^u being the spacetime metric) do not depend on time. The energy distribution is therefore 
stationary, and the total energy is 

<x> 

E = 4vr / drr W + 0' 2 + U (</>)) , (8) 
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where the prime denotes differentiation with respect to r. The field equation, 

= + (9) 

reduces to 

= A ?^M + ,v, (io) 

r d<p 

which is equivalent to 

r 

l^ + l^-u = e-2j^r. (ii) 

o 

This effectively describes a particle moving with friction in the one dimensional potential 

^(0) = lc 2 2 - (12) 

£ is an integration constant playing the role of the total 'effective energy'. It is essential 
that the potential V((p) should have the qualitative shape shown in Fig. [|, which is possible 
if the following conditions are fulfilled. First, since V"(0) < 0, it follows that u 2 should not 
be too large: 

u 2 < u£, = ^"(0) ■ (13) 

On the other hand, u 2 should not be too small, since otherwise V(4>) will be always negative. 
V{4>) will become positive for some non-zero <ft, as is shown in Fig. [|, if only 

u 2 > u 2 nm = min(2[/(0)/0 2 ) , (14) 

where the minimum is taken over all values of <fi. For the potential U(<f>) it is necessary that 

U"(0) > min(2f/(0)/0 2 ) , (15) 

since only then the set of values of uo 2 will be non-empty. The only possible renormalizable 
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interaction in the theory, U = |/i 2 |$| 2 + A|$| 4 , does not obey this condition. Thus non- 
renormalizable potentials have to be considered ||. For example, for the potential 

U(4>) = A(0 6 - a0 4 + 60 2 ) (16) 



the condition flT5p is fulfilled for any positive A, a, b, and U(<f) will have a global minimum 
at = if b > a 2 /4. We use this model potential with A = 1, a = 2 and b = 1.1 in all our 
calculations below, in which case the conditions (13), (14) require that 0.2 < u 2 < 2.2. 



III. FUNDAMENTAL Q-BALLS AND THEIR RADIAL EXCITATIONS 

If conditions (|T3] ) -([T5| ) are fulfilled, then the field equation admits globally regular solu- 
tions 0(r) with finite energy |J. The necessary condition for the energy (|8]) to be finite is 
that the potential U — > for large r, and therefore — » as r — > oo. Linearizing Eq. ( |T0"D 
around = 0, one finds that asymptotically 

= j exp {- ^"(0) - ^) r} (1 + 0(l/r)) , (17) 

where A is an integration constant. In view of ([H|) the argument of the exponent is real 
and negative, and so approaches zero exponentially fast. 

Solutions must also be regular at the origin of the coordinate system, r = 0. Since r = 
is the regular singular point of Eq. (|T0|), the solution will only be regular if it belongs to the 
'stable manifold' characterized by the local Taylor expansion in the vicinity of r = 0, 

0(r) = O + (f/'(0 o ) - o; 2 0o) r 2 + 0(r 4 ) , (18) 

where 0o is an integration constant. 

Extending the two local solutions (|T7D, (|1"8| ) to finite values of r and requiring that and 
0' for both solutions agree at some r = tq, yields two conditions for the free parameters 
A and 0o- Resolving these conditions determines a globally regular solution in the interval 
r G [0,oo). In fact, in this way an infinite discrete family of globally regular solutions 
parametrized by the number n = 0, 1, 2, ... of nodes of is obtained. The existence of these 
solutions can be illustrated by the following qualitative argument. 
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FIG. 3: The profile of <p(r) for the funda- 
mental Q-ball solution and its first two radial 
excitations. 



FIG. 4: The radial energy density v 2 Tqq for 
the fundamental Q-ball and its first two ex- 
citations. 



The parameter 0o in fllSf ) is the coordinate of the 'particle' at the initial moment of 'time', 
r = 0, when the particle velocity is zero, 0'(O) = 0. The particle therefore starts its motion 
from some point on the curve V(4>) in Fig. ^|, with the total effective energy £ being equal to 
its potential energy. Then it moves to the right, dissipating some of its energy along its way. 
For r — > oo, the particle must end up at the local maximum of the potential V (at = 0) 
with the total effective energy being zero. This can be achieved by fine-tuning the initial 
position of the particle, 0o- If V(<f>o) < 0, the effective energy of the particle will always be 
negative, and therefore it will not be able to end up in a configuration with zero energy. On 
the other hand, if 0o is such that the particle starts very close to the absolute maximum of 
V, then it will stay there for a long 'time' r, during which period the dissipation term in 
(|TT|), which is ~ 1/r, will become very small. As a result, when the particle will eventually 
start moving, its energy will be too large, and so it will 'overshoot' the position with zero 
energy. By continuity, there is a value 0o for which the total effective energy (the right hand 
side in (|Tl~D) is exactly zero for r — > oo, and so the particle will travel from 0(0) = 0o to 
0(oo) = O§. 

Next, one can fine-tune O such that the initial energy V((p ) is slightly too large, so that 
the particle first overshoots the = position, but then it hits the barrier from the other 
side, bounces back and dissipates just enough energy to finally arrive at = with zero 
energy. This will give a solution with one node of 0(r) in the interval r 6 [0, oo). Similarly 
one can obtain solutions with n > 1. To recapitulate, for each u 2 subject to (|i~3|), ([14]) there 
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is a solution to Eq. (|10"D for which 0(r) smoothly interpolates between some finite value at 
the origin and zero value at infinity, crossing zero n times in between. We shall call solutions 
with n = fundamental Q-balls. Solutions with n > are the radial Q-ball excitations. 
The family of Q-ball solutions can thus be parameterized by (uj,n). 

As the frequency uj changes in the range tt; mm < uj 2 < c<J max , the charge 

oo 

Q(u) = J drr 2 (p 2 (19) 
o 

changes from Q(a'min) = ±oo (depending on the sign of to) to Q(ciWx) = 0. This can be 
understood as follows. For uj 2 — > cu max the minimum of the effective potential V{4>) becomes 
more and more shallow and moves closer and closer to = 0. This implies that the range 
of values of the solution <f)(r) diminishes and the interval of r in which </>(r) is not constant 
decreases. Q-balls therefore 'shrink' in this limit, and the integral (|19"D tends to zero 6 . In 
the opposite limit, uj 2 — > cu min , Q-balls become large, and their charge tends to infinity. 

These considerations imply that instead of u one can choose the charge Q as the in- 
dependent parameter. Spherically symmetric Q-balls therefore comprise a two-parameter 
family labeled by (Q,n), where the charge — oo < Q < oo and the 'excitation number' 
n = 0, 1, 2, . . .. In Figs. 3-4 the profiles of the fundamental solution and its first two radial 
excitations are shown for Q = 1100. As one can see from the shape of the radial energy 
density, the n-th solution has the structure of n + 1 concentric spherical layers of energy 7 . 

IV. SPINNING Q- VORTICES. 

Our aim now is to show that Q-balls admit spinning generalizations. For this we modify 
the ansatz (|6|) according to 

$ = 0( r , tf) e ^+ iN ^ , (20) 

6 Although the exponent in ( |l7| ) decays slower and slower for ui 2 — > U"(0), the prefactor A vanishes faster 
still. 

7 To our knowledge, solutions with n > have not been reported in the literature so far. 
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where N is an integer. This produces a non-zero component of the angular momentum 
along the z-axis. Inserting this ansatz into the field equation ([5p, the problem reduces to a 
non- linear PDE for the function 0(r, Before we start solving this equation, however, it 
is instructive to consider a simpler system which effectively lives in 2 + 1 dimensions. Then 
the problem reduces to an ordinary differential equation (ODE). 

Let us pass from spherical coordinates (t, r, ip) to polar coordinates (t, p,z,ip), and then 
assume that the field does not depend on z: 



This will correspond to a vortex-type configuration, invariant under translations along the 
z-axis. The energy per unit vortex length is 

oo 

E = 2tt J p (u 2 <f> 2 + <t>' 2 + y <P 2 + U((f>^j dp . (22) 
o 

If TV ^ 0, the vortex will rotate around the z-axis, its angular momentum per unit length 
being given by 

oo 

J = J To^pdpdtp = AttcuN J pcfdp = NQ , (23) 

o 

where Q is the charge per unit length. The field equation now reads 

= # . + V-^-!^W ,24, 
p p A dcp 

For N = this reduces to (|T0|), with the replacement 2/r^l/pin the friction term. All 
arguments above still apply, hence we conclude that there exist globally regular 'Q- vortex' 
solutions and their radial excitations with finite energy per unit length. These solutions 
display the behaviors qualitatively similar to those shown in Fig. || 

Let us now consider solutions with N > 0. As can be seen from (|22D , the energy for such 
solutions will be finite if only 0(0) = 0, such that the boundary condition for small p is now 
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different. The power series solution to (|24j) in the vicinity of p = reads 



Bp N + 0{p N+1 ) 



(25) 



where B is an integration constant. The asymptotic behavior for large p is 

A 



VP 



cxp 



■y/(U»(0)-U*)p\(l + O(l/p)) . 



(26) 



One can use essentially the same qualitative considerations as in the preceding section to 
argue that globally regular solutions to ( p4|) with such boundary conditions exist, provided 
that uj still fulfills conditions (|l"3"D, fl!4]) . These solutions can be obtained by numerically 
extending the asymptotics (pop , (p6|) to finite values of p and adjusting the free parameters 
A, B to fulfill the matching conditions at some p = p . 
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FIG. 5: The spinning iV = 1 fundamental 
Q-vortex solution and its first two radial ex- 
citations; Q = 150. 
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FIG. 6: The first three spinning excitations 
of the fundamental Q-vortex solution; Q = 
60. 



The conclusion is that there exists a family of globally regular spinning Q-vortex solutions. 
These solutions can be labeled by (Q,n,N), where Q is the charge per unit vortex length, 
n = 0,1,2,... is the radial 'quantum' number (the number of nodes of <f)(p)) and iV = 
0,1,2,... is the rotational 'quantum' number 8 . For fixed Q, the energy per unit length, 
E(Q,n, N), depends on both n and N, while the angular momentum is determined only 
by the value of iV as J = NQ. The profiles of <fi(p) for several low-lying excitations of the 



Spinning Q- vortices with n = have also been found in Ref. (H 
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fundamental Q-vortex are shown in Figs. |5|-|| 



V. SPINNING Q-BALLS 

Having considered the simpler problem in 2 + 1 dimensions, we now return to our main 
task of finding rotating Q-ball solutions in 3 + 1 dimensions. Although, these two cases are 
qualitatively somewhat similar, the 3 + 1 dimensional problem is technically more involved, 
since it requires solving a non- linear PDE. With $ = <j>(r,$) e iuJt+iNip the field equation 
reduces to 

— + -— + — — + cos ^ 9 _ N<2 +J*\(h - dU ^ (27) 
The energy, E = J Tood 3 x, reads 

°° n 2 2 

E = 2vr [drr 2 [ dd sintf f u 2 ^ + (<9 r 0) 2 + \{d^<f)) 2 + f . t - + U(<f>) J , (28) 
J J \ r z r l sin v J 

o o 

the charge 

OO 7T 

Q = 2uj J <p 2 d 3 x = 4nu J drr 2 J d$ sin$0 2 , (29) 

o o 

and the angular momentum 

J = Jt 0ip d 3 x = NQ . (30) 
Finiteness of the energy requires that 

as r-^0oroo. (31) 

The asymptotic behavior of the solutions in these limits can be easily determined, since for 
small one has dU/d(p ~ U"(0)4>, such that equation fl27|) actually becomes linear. The 
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variables then separate, implying that the most general asymptotic solution has the form 

oo 

0M) = 5^/,(r)if(oo8tf). (32) 

l=N 

Here, Pf* (cos $) are the associated Legendre functions. At the origin, the radial amplitudes 
fi(r) are 

/,(r) = Qr l + 0(r l+1 ) , (33) 

while at infinity 

fi(r) = y exp{-v/(f/"(0)-^)r} (1 + 0(l/r)) , (34) 

C\ and Ai denoting integration constants. 

Eqs. (j32^)-(|34]) have been obtained by linearizing the field equation in the vicinity of r = 
and r = oo, in which case modes with different values of the quantum number I decouple. 
Our strategy to construct solutions in the whole space is to employ again the partial wave 
decomposition (|32]). This is always possible, since the associated Legendre functions form a 
complete set. However, since for arbitrary values of r the equation is non-linear, harmonics 
with different values of / will no longer decouple. 

Since Eq. fl2~T|) is symmetric with respect to reflections in the xy-plsme, — ► 7r — it 
follows that if 0(r, i?) is a solution, so is 0(r, ir — In addition, — (f>(r, $) is also a solution, 
since the field equation contains only odd powers of <p. The associated Legendre functions 
P/^cost?) are even/odd with respect to i? — > 7t — i9 for even/odd values of l + N, respectively. 
As a result, half of the terms in the mode expansion (|32|) will change sign under the reflection, 
while the other half will stay invariant. Since 0(r, tc — $) must also be a solution, it follows 
that either all odd or all even terms in the mode expansion must vanish in order to have 
either <ft(r, tc — -&) — cf)(r,d) or <f)(r,7t — i?) = — <p(r, The conclusion is that for a given 
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value of N there are two solutions with different parities P: 

oo 

P = +1 : 0(r,tf) = ]T/ fc (r)P^ +2fe (costf) , (35) 

k=0 

oo 

P = -1 ; 0( r ,tf) = ^/ fc (r)P^ +2fc+1 (cos^) . (36) 



k=0 



With our choice of the potential (jI6|), the field equation fl2T|) contains cubic and quintic 
non-linearities. In view of the completeness of the associated Legendre functions, their prod- 
ucts can be expressed in terms of their linear combinations, for example, (Yl fi^t 1 i x )) 5 = 
^ / 0'jPf(x), with the coefficients aj determined by the //'s. As a result, inserting fl3"5]), ( |3"6] ) 
into (|27|) we obtain 

oo 

J2v k [fs(r)]P% +k (cos$) = 0. (37) 

k 

Here T>k[f s { r )] are second order differential operators acting on the radial amplitudes / s (r), 
and the sum is taken over all odd/even positive k for odd/even solutions, respectively. This 
equation is equivalent to the infinite set of ODEs 

T>k[fs(r)} = . (38) 

Now, we truncate this system by setting all amplitudes / s with s larger than some maximal 
value Z max — N to zero and by discarding all equations with k > Z max — N. The indices in 
(f35|) then vary only in the finite limits 

k,s = 0, 1, ... , Z max - N . (39) 

As a result, we end up with a finite system of ODEs. This procedure is sometimes called 
Galerkin's projection method. It is natural to expect that if Z max is large enough, the resulting 
approximate solutions will be close to the exact ones. To illustrate that this is indeed the 
case, we show in Table I the energy and charge of the solution of the truncated system with 
uj 2 = N = P = 1 for several values of Z max . As one can see, the energy and charge indeed 
approach some limiting values with growing / max . It actually seems to be sufficient to take 



16 



TABLE I: Parameters of the solitons versus the truncation parameter I 



^max 


1 


3 


5 


7 


9 


E 


230.70 


217.16 


207.74 


205.58 


205.36 


Q 


211.46 


200.68 


189.43 


186.99 


186.73 



into account only the first 3-5 lowest harmonics in order to get a reasonable approximation. 

In the next step we construct solutions for a fixed charge Q in different N sectors. More 
precisely, the numerical solutions to Eqs. ( |3~8"D have been obtained with Matlab®'s ODE 
solvers by utilizing the shooting method. The asymptotic solutions (|33|), (|34]) were used to 
start the integration at r = 0.01 and at r = 10 toward the matching point, whose position 
has been varied between r = 3 and r = 6. The matching conditions imposed at this point 
determine the values of the constants C\ and Ai in fl3"3] ) and fl3"4|). The typical matching 
error was found to be less then 10~ 16 . The profiles of even and odd rotating solutions with 
N — 1 are shown in Figs. 0-[l^. As one can see from these plots, the distribution of the 
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0" 00 



FIG. 7: 0(r, i?) for N = 1, P = 1. FIG. 8: T 00 for N = 1, P = 1. 

energy density T 00 is strongly non-spherical. It has the structure of a deformed ellipsoid for 
P = 1, and that of dumbbells oriented along the rotation axis for P = —1; in the latter case 
the energy density vanishes in the equatorial plane. The energies of the first three rotational 
excitations of the fundamental Q-ball are given in Table |I[ As one can see, the energy of the 
first excitation exceeds the ground state energy by about 20%, the next excitation lying again 
about 20% above. This is in agreement with the expected properties of spinning excitations 
of a single soliton. In summary, spinning Q-balls comprise a two-parameter family labeled 
by the values of (Q,N). We have also found evidence for the existence of spinning radial 
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TABLE II: Parameters of the rotating solutions with Q = 410. 



t\tP 

Jy 


7 

'max 


hi 




U 


U 


qa7 on 


u. ( ouyy 


1 + 


9 


378.36 


0.86017 


1- 


8 


442.24 


0.99164 


2+ 


10 


433.94 


0.96258 


2" 


9 


505.66 


1.13284 


3+ 


9 


473.59 


1.05277 


3" 


12 


528.30 


1.31037 



excitations with n > 0. However, the construction of such solutions is somewhat more 
involved and we refrain from presenting them in this paper. 




" 

FIG. 9: <p(r, &) for the N = 1, P = -1. FIG. 10: T 00 for the N = 1, P = -1. 

VI. CONCLUSIONS 

The aim of this paper was to find an example of spinning solitons in Minkowski space. 
We have considered the model of a complex scalar field with a non-renormalizable self- 
interaction. In the spherically symmetric sector this model contains non-topological solitons, 
the Q-balls. In addition, we have found an infinite discrete family of radial Q-ball excitations, 
parameterized by the number of nodes n of the radial field amplitude. Such excited solutions 
have not been reported in the literature before. 

In a second step we have analyzed cylindrically symmetric solutions with explicit har- 
monic dependence on the azimuthal angle, exp(iNip), which we call spinning Q- vortices. For 
such solutions there is a non-zero component of the angular momentum along the z-axis, 



18 



J = NQ, where Q is the charge per unit vortex length. In addition, these solutions ex- 
hibit radial excitations parameterized by an integer n. As a result, such spinning solutions 
comprise a three-parameter family labeled by (Q,N,n). 

Finally, we have considered the full 3 + 1 dimensional problem. We have used a version 
of the spectral method by expanding the field with respect to the complete set of associated 
Legendre functions. We reduced the PDE to an infinite system of radial ODEs, and then 
truncating this system at finite order. The parameters of the solutions for the truncated 
system converge rapidly to some limiting values as the truncation parameter grows. By 
keeping the charge Q fixed, we have obtained the lowest rotational excitations of the funda- 
mental Q-ball solutions. The angular momentum of these solutions is quantized, J = NQ, 
the energy increases (but not very rapidly) with the angular momentum. As a result, these 
solutions can be viewed as describing spinning excitations in the one-soliton sector rather 
than orbital motion in a many-soliton system 9 . To our knowledge, these spinning Q-balls 
provide the first explicit example of spinning solitons in Minkowski space in 3 + 1 dimensions. 
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